function [g0dd g0UU g0det] = h2basemetric(params)

  % the computed quantities have the asymptotic behaviour pulled out as follows:

  % g0dd -- 1/z^2
  % g0UU -- z^2
  % g0det -- 1

  T = params.TempOffset;
  z = params.Cz;


  g0dd  = zeros(params.Nz, 4, 4);
  g0UU  = zeros(params.Nz, 4, 4);
  g0det = zeros(params.Nz);

  g0dd(:,1,1) = (-1)+T.^3.*z.^3;
  g0dd(:,1,4) = -1;
  g0dd(:,4,1) = -1;
  g0dd(:,2,2) = 1;
  g0dd(:,3,3) = 1;
  

  g0UU(:, 1,4) = -1;
  g0UU(:, 4,1) = -1;
  g0UU(:, 2,2) = 1;
  g0UU(:, 3,3) = 1;
  g0UU(:, 4,4) = 1+(-1).*T.^3.*z.^3;

  g0det = -1./z.^8; 
end
